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Abstract 

For the mathematical modeling of highly non-eq\iilibrium and non- 
linear processes in a tornado in this paper a new approach based on 
nonlinear equations of momentum transfer with function of sources and 
sinks is suggested. In constructing the model thermodynamic description 
is used, which is not entered before and allows discovering new principles 
of self-organization in a tornado. This approach gives fairly consistent 
physical results. This is an attempt to answer some fundamental ques- 
tions concerning the existence of a tornado based on the created model 
and numerical results. 
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1 Introduction 

Tornado (Spanish Tornado "Whirlwind"), twister - rapidly rotating at- 
mospheric vortex that occurs in cumulonimbus (thunderstorm) cloud and 
spreads down to the ground, in the form of cloud arms or trunk of tens or 
hundreds of meters. 

This phenomenon exists for millions of years and it is also observed on 
other planets with atmosphere. But until now it is a little studied and hard 
to be measured owing to the presence of large air velocities near the funnel 
(kernel, trunk) and, consequently, owing to the inability of the measuring 
equipment to withstand powerful external influences [2j. 

From a purely dynamical point of view, and tornadoes occur due to the 
strengthening of existing localized vortices by an external flow of momentum [4] 
Supercell thunderstorm models [H[3l[6t [TTtll3tll7j presented in the last forty 
years, have many similarities in their basic structure. However, in contrast 
to the purely dynamic approach to work, beginning with article [TT] all mod- 
els have intensive upstream, which cross the vertical wind shear, increases 
its speed and swirl with height. Thus, these models are the convective and 
they are based on the model of mesocyclone, which formation has been 
just recognized as the most useful fact to supercell identification [4] . At the 
same time, they differ in the detail level from the occurring processes and 
as a consequence in the initial and boundary conditions that are necessary 
for the thunderstorm supercell occurrence. Numerical models of supercell 
thunderstorms [HI I19j are based on the dominant mesocyclone. The hydro- 
dynamic model [T9], which takes into account a set of cross-produced pair 
and other interactions and forces, is particularly interesting. 

Observations of the tornado have a rich history, provided by many papers 
only for the 20th century. Brooks was the first observer, who put forward 
generally accepted assumption, that the funnel is a part of the parent cloud, 
the structure and dynamics of which represent a small tropical storm and 
having a helical structure [2]. Numerous observations of the parent cloud 
indicated the presence of long vortices in the horizontal plane in them; the 
vertical poles (funnels) are the continuation of which [18] . This fact has no 
explanation. In 1951 in Texas during a tornado the funnel, passing over the 
observer, rose, and its edge was at the height of 6 meters with the inner 
cavity diameter of 130 meters. The wall thickness was the size of 3 meters. 
Vacuum in the cavity was absent, because it was easy to breathe during 
its passage. The walls was extremely fast spinning (Justice, 1930). In the 
monograph Flora [5j notes that "the distinction between the strongest winds 
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in the body of the funnel and the stationary air on its periphery is so sharp, 
that it causes a number of damaging effects." 

Observations of the actual tornado, therefore, indicate a strong non- 
linearity and non-equilibrium of processes in atmosphere during the for- 
mation and existence of a tornado, that does not allow to create the perfect 
model of this exotic phenomenon. 

In the framework of the study of this unusual natural phenomenon the 
following questions are need to be answered: 

1. Under what conditions in the atmosphere the appearance of a tornado 
happens? 

2. What causes the existence of distinct lateral boundaries of the tornado? 
Why don't these boundaries spread in time? 

3. What are the conditions for the existence of dissipative structures 
observed in the tornado - a set of organized vortices? And what are their 
conditions of decay? 

4. What causes the boundedness of tornadoes in height? 

5. What are the conditions for the existence and stability of stationary 
and others possible modes. 

6. What determines the appearance of the tornado core (trunk) with 
significantly higher velocities? 

For the mathematical modeling of highly non-equilibrium and nonlinear 
processes in a tornado authors propose the approach based on the nonlinear 
equations of momentum transfer with the model sources and sinks function. 
This approach can be assigned to the problems with peaking considered by 
Academician A. A. Samarskii [51 114t \15\ I16j . For the first time the thermo- 
dynamic description was used to identify new principles of self-organization 
in the atmosphere in the model specification not entered before. 

2 Model 

On the basis of the equations of momentum with a nonlinear function of 
sources and sinks the nonlinear hydrodynamic model of a strongly nonequi- 
librium processes in the atmosphere during an intense vortex formation was 
constructed by the authors. A thin layer of a unit volume, parallel to the 
surface of the earth, is allocated in the tornado. In this layer there are 
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sources of the movement and its sinks depending on the horizontal velocity 
vector and its modulus. The problem is solved under the assumption of 
isothermality of the layer and its vicinity. Temperature and other air char- 
acteristics changes with the altitude. "Layered" model is used to simplify 
the modeling of vortices existing in the three-dimensional region of space, 
as well as to facilitate numerical calculations. 

Non-linear momentum source in the atmosphere in some cases leads to 
the regime with peaking - sharp increases in the velocity amplitudes during 
some short-time. The development of this mode leads to a self-organization. 
In this formulation, the formation of dissipative structures localized in a 
space, is possible. In these structures the speed may increase limited (or 
unlimited) . 

Authors using this model based on the thermodynamic approach, attempt 
to describe the observed physical phenomena and to explain the mechanisms 
of nonlinear layered momentum transfer in a tornado in case of the nonpo- 
tential flow. 

In the non-equilibrium thermodynamics it is accepted to characterize the 
processes within the system under the influence of the external environment 
by the so-called entropy production a* per unit volume of the layer. There 
are also other local thermodynamic characteristics - the external flow of 

entropy a'^ and the rate of change of entropy S, equal to their sum. It 
is believed that in the self-organization systems full change of the entropy 

decreases with time: 5* < 0. 

This approach allows one to record the system of equations for the velocity 
components in the case of an incompressible fluid in dimensionless form as 
two-dimensional Kuramoto-Tsuzuki equation[10j for the atmospheric layer: 

(1) — = u*, (1 + ^Cl) + ^ j + q*<!> - a*, (1 + ^C2) 

where ^ = '&% + iiD*; Ai = 1^2/1^1 related to the viscosity, A2 = a\la\ due to 
sinks. The superscript " *" means the dimensionless form of the parameters 
on certain scales, identified at the observing of a tornado. 

3 Results 



The results of numerical simulations give a good resemblance to the observed 
physical phenomena. Approach allows us to explain nonlinear mechanisms 
layered momentum transport in the atmosphere, and also origin, evolution 
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and decay of large, medium and small atmospheric vortices as dissipative 
thermodynamic structures. Competition of pulse increment and propagation 
processes in viscous medium leads to the appearance of a linear size - the 
spatial diameter of self-organized structure - Iq on the border of which the 
rate of change of entropy changes sign. You can set the dependence of the 
initial conditions from humidity, which varies with altitude. 

Humidity, which determines the radial component of the emerging spiral 
wave, in the lower layers at its reduction leads to a significant depletion of 
the vortex structure around the trunk and even their disappearance, and, 
thus, to more realistic result when the trunk is only visible between the 
earth and cloud formed the tornado (Fig. [T]). 




Figure 1. The entropy production in the tornado basin in 
layers at different heights from the ground up to 3.5 km. 

The Table [T] shows the spatial characteristics of the comparatively small 
size tornado defined by numerical methods for two variants of the approx- 
imation of turbulent flows: nonpotential (NP) and potential (P). Nonpo- 
tential flow means that the function of the pressure gradient relates to the 
velocity in terms of the sources and sinks. Determination was made on the 
velocity (i?*) and pressure gradient (Vp*). The obtained characteristics are 
in a good agreement with observations of real tornado, that is according to 
the authors confirm the conformity between theory and experiment. 
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Thus, the thermodynamic model is fundamentally different from the con- 
vective models, and allows you to construct model that reflects some of the 
principal features of the real observed tornado. This is such features as the 
presence of the trunk almost perfectly circular in cross section, the existence 
of significant gradients of velocity and pressure on the edge of the tornado 
basin and in the area of the vertical vortex tube [5], along with a complete 
absence of rotation and discharge in the center of the tornado [7J. The con- 
sidered approach allows us to formulate conditions of self-organization in 
a tornado established by the authors in present. This approach allows us 
to expand the problem in the asymptotic to conditions of formation of cy- 
clones in the atmosphere. Also some physical restrictions for the horizontal 
velocity, the intensity of the sources and sinks and the stability condition of 
dissipative structures are specified. These restrictions can not be obtained 
explicitly using the convective and other models. 

Such a nonlinear hydrodynamic model based on a thermodynamic ap- 
proach can be used to describe cyclones, tornadoes, and other exotic natural 
phenomena. 

3.1 Condition for the stability of vortex flows in a tornado. 

Stability analysis of solutions of the Kuramoto-Tsuzuki equation [10] yields 
the following result: 

(c? + 1) A;^ + 2 (1 + C1C2) e>0,k = tt/Iq. 

In order to satisfy the inequality for any value of k it is necessary to limit 
constants ci,C2: — 1 < C1C2 < 1. This condition is the criterion for the 
stability of turbulent structures in atmospheric vortices. Therefore, it is a 
condition of release of such vortices in the steady state. 



Table 1. Spatial characteristics of tornado with a diameter 
pool Iq = 500 m, \m\ = 1, the comparative table. 



Flow type 


NP 


NP 


P 


Characteristics 


1?* 


Vp* 


Vp* 


Self-organization zone diameter d, m 


439 


459 


426 


Inner core diameter, m 


11 


9 


9 


Outer core diameter, m 


20 


20 


21 


Width of the pressure equalization 








ring to atmospheric pressure, m 


30 


20 


37 
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The emergence of a stationary core in a plane layer of the tornado - 
space localized figure - is the result of sharpening and confirmed by the 
work[12j, which describes the unusual properties of the presence of large 
pressure gradients and velocities in the tornado. 

Real tornado in this case goes to steady state and has the core (trunk) 
formed in a planar layer, spatial characteristics of which are shown in Ta- 



As an intermediate generalization of this work we present currently in- 
stalled thermodynamic and hydrodynamic conditions of self-organization. 



3.2 The thermodynamic conditions of self-organization. 

The sum of the external (a^) and internal thermodynamic fluxes (cr*), is 

characterized by the full rate of change of entropy S- The decrease in en- 
tropy with time corresponding to processes of self-organization of vortex 
structures. The transition from the reduction of the entropy to its growth 
allows us to define the boundary of the vortex tornado basin, in which stable 
self-organized vortex systems exist. As follows from Fig. [2j the entropy of 

each tornado vortex layer decreases with time: S < 0. This is achieved with 
a value of sources of momentum - the reversible flows of entropy < 0, 
which are compared with the entropy production \a^\ > a*. The condition 

of the limited lateral tornado is a value S = 0, that is cr'^ -|- cr* =0 on the 



vortex- formation boundary (Fig. ^). Outside the vortex tornado basin en- 




tropy increases (5 > 0). The zone of self-organization decreases over time, 
which finally leads to the collapse of a tornado. 



3.3 Hydrodynamic conditions of self-organization in a tornado, 

established for the first time, follow from the thermodynamic conditions. 
Boundedness of tornado in height is owing to the absence of vortex forma- 
tion due to the high viscosity. The transverse size of the tornado is caused 
by the competition between the rate of a velocity and a momentum distri- 
bution in terms of the viscosity. Positive feedback between the projections 
of horizontal velocities leads to sharpening. The intensity of the sources and 
sinks of movement causes the limitation on the modulus of the horizontal 
velocity: > 4g*/a^(l — 4c|/9). Numerical calculations of such problem 
with the sharpening have shown the possible existence of enormous speeds 



ble 1. 




BYSTRAI G.P., LYKOV LA. 



7 




Figure 2. The total change in entropy in a plane layer of 
the tornado basin (a), its projection in the form of isolines 
(b) for the central region of the tornado that characterize 
the region of self-organization (c) and the boundary of the 
localized vortex structures for the calculation of the total 
area. 



of air movement directed to the nucleus. Also the approach allows us to cal- 
culate the spatial distribution of pressure in the tornado basin with change 
of height by the comparison with external pressure. 

Restriction on the constant of the Kuramoto-Tsuzuki equation for the 
emergence of self-organization is a condition: A2 < 02'^ /a]^. The left twist- 
ing condition in the topological charge is: m = — 1, right - m = +1. Re- 
striction on the modulus of the topological charge for the formation of a 
tornado: |m| < 2. The transition to the values of the topological charge 
|m| > 2 firstly, leads to the huge cyclone model, i.e. atmospheric structures 
without a trunk with a significantly lower velocities and, secondly, to the 
emergence of "eye" - the cyclone core. 

The localized structure described and obtained in the framework of the 
sharpening problems on the basis of the thermodynamic approach is not 
unique. Similar effects occur in the thermal and diffusion processes, ferro- 
magnetic materials. But in this article, self-organization in a regime with 
peaking for hydrodynamic problems was described by authors for the first 
time. This approach provides sufficiently uncontroversial physical results. 
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Thus, in this paper, the authors on the basis of the created model and 
the results of numerical calculations attempt to answer on all posed in the 
beginning of this article. 
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